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We estimate the amount of the missing baryons detected by the Planck measurements of the 
cosmic microwave background in the direction of Central Galaxies (CGs) identified in the Sloan 
galaxy survey. The peculiar motion of the gas inside and around the GGs unveils values of the 
Thomson optical depth tt in the range 0.2-2 x 10“"^, indicating that the regions probed around 
CGs contain roughly half of the total amount of baryons in the Universe at the epoch where the 
CGs are found. If baryons follow dark matter, the measured rxs are compatible with the detection 
all the baryons existing inside and around the GGs. 

PACS numbers: 98.52.Eh, 98.62.Py, 98.70.Vc, 98.80.Es 


Introduction. The interplay between baryons and dark 
matter is a key problem in cosmology and galaxy for¬ 
mation. Understanding the distribution of baryonic and 
dark matter in galaxies, groups, and clusters of galaxies 
is an essential step towards the full picture of how these 
objects form and evolve. It is well known [H-Q that only 
about 10% of all baryons in the Universe reside in the 
form of stellar mass, while other 90% resides in a diffuse, 
mostly undetected component, if it were not for recent 
UV spectroscopic measurements of un-collapseR diffuse 
gas in the direction of certain QSOs, [see, e.g.,y, y, and 
references therein.] Nowadays there is an ongoing debate 
[id on whether a significant fraction of the latter com¬ 
ponent is present in the circumgalactic medium around 
halos or if instead most of the gas has been expelled 
or never accreted due to feedback processes like galactic 
winds or Active Galactic Nuclei activity. This connects 
the missing haryon issue with the complex problems of 
feedback and galaxy formation. Certainly a more com¬ 
plete census of the baryon distribution in the universe 
would be of great relevance in this context. 


The kinematic Sunyaev-Zeldovich effect [hereafter kSZ, 
[I^ describes the Doppler shift in frequency induced 
on photons of the cosmic microwave background (here¬ 
after CMB) photons after they scatter off free electrons. 
This Thomson scattering induces frequency independent 


brightness temperature anisotropies in the CMB, which 
are given by |13j [ 


'5TkSz(n) = -To J dl aTUe • n^ ~ -Tq tt 

( 1 ) 

In this expression, the integral tt = ctt / dl is con¬ 
ducted along the line of sight (LOS) given by h. We 
have made the approximation that the typical correla¬ 
tion length of LOS velocities (given by v ■ h) is much 
larger than the density correlation length, in such a way 
that the LOS velocity term may be pulled out of the kSZ 
integral. This is justified by the results of Planck Collab¬ 
oration Int. XXXVII IJ], who find a typical correlation 
length of peculiar velocities of 20-40 h~^ Mpe, well above 
the typical galaxy correlation length (^ hh~^ Mpc). The 
expression above also shows that the kSZ constitutes an 
integral over the electron momentum, independently of 
the temperature, and thus counts all electrons in the bulk 
flows, regardless they belong to collapsed structures or 
not. 

Large scale, bulk matter flows were detected via the 
kSZ effect firstly by 151 and more recently by [l^. We 
build upon the latter work and extract physical con¬ 
straints on the amount of baryons contributing to the 
kSZ signal and the implications of those measurements 
in the problem of the missing baryons. As in 0 , we 
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FIG. 1. Left panel: Measured pairwise peculiar momenta 
under different aperture sizes for the SEVEM map. Right 
panel: Measured (r) correlation function under dif¬ 

ferent aperture radii (6iAp) for the SEVEM clean map. In both 
panels the green solid lines provide the expectation from the 
GALAXY mock catalogue, while the dash-dotted lines corre¬ 
spond to expectations found after averaging over 100 Gaus¬ 
sian simulations. In both cases, we assume tt = 5 x 10~®. 
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FIG. 2. Left panel: Probability distribution for the optical 
depth Tt inferred from the (r) correlation function 

(solid lines) and the kSZ peculiar momenta (dashed lines), 
under different apertures, for the SEVEM map. These estimates 
of Tt are obtained after fitting the data to the prediction in¬ 
ferred from the average of the Gaussian simulations. Right 
panel: Same as left panel, but using the GALAXY N-body pre¬ 
diction as a reference, rather than the Gaussian simulations. 


use the latest Planck data release (pR2) available at the 
Planck's Legacy Archive server [16|, the Central Galaxy 
Catalogue (hereafter CGC) obtained from the Sloan Dig¬ 
ital Sky Survey SDSS/DR7 0, and a mock catalogue 
of central galaxies obtained from the Millennium numer¬ 
ical simulation [l^, to which we henceforth refer as the 
GALAXY mock catalogue. In this work we quote results for 
the Planck SEVEM CMB map, although similar results are 
obtained for all the other cleaning algorithms. Our esti¬ 
mates of the kSZ temperature of Eq. (ED are taken from 
M and are obtained after applying aperture photome¬ 
try (hereafter AP) in the direction of the CGC members. 
We defer the reader to that paper for a more detailed 
description of the data under use, and a deeper discus¬ 
sion on the measurements of the kSZ, while in this work 
we concentrate upon the physical interpretation of such 
measurements. In particular, we set constraints on tt 
and the missing baryons from our AP measurements of 
JTksz in Eq. [TJ 

Methodology and results. The two statistics presented 
in [l^ yielding kSZ evidence are the kSZ pairwise pe¬ 
culiar momentum (hereafter pksz(?')) and the cross cor¬ 
relation function of the kSZ temperature and the recov¬ 
ered radial peculiar velocity (hereafter (r)). The 

Pksz (f) reflects the gravitational infall in pairs of galaxies 
and is computed from the sum 


Pkszir) = - 






( 2 ) 


i<j id 


where the weights Cij are given by [If 


(d - g)( 1 + cos6<) 

2 +r^ - “^riTj cos t 


( 3 ) 


Following the convention of [l5j, in this equation and 
Tj are the vectors pointing to the positions of the i-th and 
j-th galaxies on the celestial sphere, and Vj are the co¬ 


moving distances to those objects, and 
refers to the distance vector for this pair of galaxies. 
The hat symbol (f) denotes a unit vector in the direc¬ 
tion of r, and 9 is the angle separating and i-j. Note 
as well that fij refers to the direction of the difference 
vector, i.e. (g — rj)/\ri — r^j. On the other hand, the 
w^Tv function is computed as the spatial correla¬ 
tion function between the measured kSZ temperature 
estimates and the linear theory predictions of the ra¬ 
dial peculiar velocity of the CGs (?i‘'®‘^), w^’^^"°(r) = 

WjUij , where as before the i,j 
sub-indices refer to CGs and Wi, Wj to their weights. For 
both statistics, the STjS correspond to the AP estimates 
of the kSZ fluctuations of Eq. ED- We refer again to [lj| 
for details on how linear theory predictions of GG radial 
peculiar velocities are obtained. Figure [T] displays the 
measured Pksz{i") and (r) functions under 5, 8, 

and 12 arcmin aperture radii for the SEVEM clean map 
of Planck (very similar estimates are found for other 
GMB clean maps). In order to estimate the amount of 
free electrons giving rise to these signals, we compare 
these measurements with predictions obtained from (i) 
the GALAXY mock catalogue of central galaxies, and (ii) 
a suite of 100 Gaussian simulations of the matter den- 
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sity contrast field that is then inverted into a peculiar 
velocity field by means of the linearized continuity equa¬ 
tion. The Gaussian simulations of the density contrast 
are generated from a dark matter linear power spectrum 
compatible with Planck's cosmology, and have no power 
at wavenumbers k > 0.15/iMpc“^ since those are re¬ 
garded as non-linear scales. The first approach provides 
velocities for the halos hosting the CGs, while the sec¬ 
ond computes an estimate of the smooth, linear peculiar 
velocity in a region surrounding each GG. Following the 
approximation in Eq. o, green curves in Fig. [1] provide 
predictions for Pksz(?') and (r) from the GALAXY 

mock (solid lines) and the Gaussian simulations (dash- 
dotted lines), with a choice of tt = 5 x 10“® for display 
purposes. We hnd that, while for (r) the two 

predictions just differ in a ^ 15 % amplitude factor, for 
Pksz(?') they differ both in shape and amplitude. Since 
the kSZ is built upon all electrons present in the volume 
sampled by the aperture photometry, and not necessar¬ 
ily bound to the CGs, we expect the real signal to lie 
between the two predictions. 

Following the template fit approach used in 0 that 
conducts minimum hts of the data to the predictions 
while assuming Gaussian and correlated errors, we ob¬ 
tain the density probabilities for tt as given in Fig. [2] 
When fitting to the predictions inferred from the GALAXY 
catalogue, we obtain tt estimates frompkSz(?') measure¬ 
ments falling in the (0.1-1.1) x 10“"^ range, in slight 
tension (between 1 and 2 a low) compared to estimates 
from (r). A better consistency among tt esti¬ 

mates is found when fitting to the predictions provided 
by the Gaussian simulations (see left panel of Fig. [ 5 ]). 
The Pksz{r) and (r) measurements should yield 

similar tt estimates, and hence Fig. suggests that the 
motion of electrons in and around CGs is better described 
by the Gaussian simulations, which we adopt in subse¬ 
quent analyses. Note that this choice is not crucial for 
our discussion below. 

We next estimate the baryon fraction giving rise to the 
observed tt amplitudes under different apertures. For 
any line of sight throughout any CG, all electrons in¬ 
side a cylinder whose base radius is given by the angular 
aperture (0ap) should contribute to the signal. The LOS 
depth of this cylinder corresponds to the comoving length 
of typical velocity correlation (Lcorr), and is suggested 
by Fig. [T] to he in the 20-40Mpc range (comoving 
units). We find however that the baryon fraction should 
not depend upon Lcorr for any given aperture: 


{fh)z{dAp) 


TT{9AP,z)/{(TTne{z)L 

corr /(1+^)) 

l//vol(2) 



^ {tt{9ap,z))z (ncGiz)dl(z))^, ( 4 ) 

where fvoiiz) — ncciz) Tr{r{z) 9 ap)^Lcoii: is the fraction 
of the comoving volume sampled by cylinders centred 
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FIG. 3. Top panel: Baryon fraction versus aperture an¬ 
gle, as inferred from (r) (filled circles and squares) 

and Pksz(r) (filled and empty triangles) measurements. Red 
(green) color denotes positive (negative) values. Bottom 
panel: Recovered profile of gas overdensity from (r) 

measurements and for Lcorr = 38A“^Mpc (filled symbols). 
In both panels, we compare our data to the predictions for 
dark matter [as provided by [13, which are given by the blue 
lines. 


upon GGs, Lcorr is the LOS depth of the cylinder in 
comoving units, he( 2 ) is the average physical electron 
number density at redshift z (whereas the “0” subscript 
denotes at present), ncciz) is the GG comoving number 
density, and r(z) and dA{z) are the comoving and angular 
distances to redshift z, respectively. The brackets 
denote redshift averages, and thus our previous measure¬ 
ments of Tt correspond to rT(0Ap) = (tt(^ap, z))z- Note 
that Eq. O holds as long as there is no overlap be¬ 
tween cylinders along two different LOSs: we have veri- 
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fied that such overlap is negligible for apertures smaller 
than 13 arcmin, (it affects less than 1 % of pairs falling 
in the first distance bin at ^ 9/i“^Mpc, and this ratio 
is still smaller for the other distance bins). We remark 
that this definition of baryon fraction computes the ratio 
of detected baryons around CGs to the total amount of 
baryons. Eq.|4] shows that this baryon fraction depends 
on the CGC completeness via the product of the average 
GG number density times the average tt measured in 
and around those galaxies. 

The top panel of Fig. [3] shows that the baryon fraction 
increases with the aperture, amounting to 45-55 % of the 
total of all baryons for 0 ap ~ 8-10 arcmin. For larger 
apertures, our measurements become compatible with 
noise and oscillate around zero. The blue lines provide 
the prediction for a scenario where the baryons trace per¬ 
fectly the distribution of dark matter. These predictions 
follow the approach of [2^ , which compute the dark mat¬ 
ter - central halo correlation function ^h,m(?')(oc Spni{r)) 
which is dependent of both mass and redshift of the ha¬ 
los, and requires some linear bias relation b{M, z) that 
we adopt from 2l|. These predictions compnte, as in 
Fq. (H, the ratio between the mean dark matter over¬ 
density inside the cosmological volume sampled by the 
cylinders centred upon the GGs and the inverse of the 
fraction of the total cosmological volume sampled by the 
same cylinders. We thus need to integrate ^h,m(?’) inside 
the cylinders, and as shown in the plot these predictions 
are independent of the choice of Tcorr- We can also incor¬ 
porate the impact of our AP filter in our prediction, after 
taking into account the redshift and mass distribution of 
our GG sample for each aperture ^ap- We find that our 
measurements of the baryon fraction are compatible with 
this prediction, although slightly falling on the high am¬ 
plitude side: results from (r) are given by filled 

red circles (filled green squares) when positive (negative), 
while filled red triangles (empty green triangles) corre¬ 
spond to positive (negative) /b estimates inferred from 
the Pksz {r) . A template fit to the dark matter prediction 
(given by the blue lines and performed as in Eqs. 11-12 
in 0, thus accounting for covariance among different 
apertures) yields an amplitude of Adm = 1-15 ± 0.30 for 
baryon fraction estimates obtained from the (r), 

and Adm = 1.20 ±0.22 frompksz(c) measurements of the 
baryon fraction. In both cases we compare these observ¬ 
ables to the prediction from Gaussian simulations. As 
mentioned above, tension at the 1.6 cr level arises when 
comparing fits of the two (r), pksz(c) observables 

to the N-body GALAXY prediction: Adm = 1-49 ± 0.37 
vs Adm = 0.85 ± 0.15 from (r) and Pksz{r) data, 

respectively. We remark that calibration errors in GMB 
maps would impact these amplitudes while leaving the 
significance of the kSZ signal unchanged. 

Alternatively, in the bottom panel of Fig. [3] we com¬ 
pare the gas overdensity values inferred from the rx mea¬ 
surements from our " (r) observations and Leon = 


38/i“^Mpc (filled symbols) with the predictions of dark 
matter according to (blue dot-dashed line), for dif¬ 
ferent choices Tcorr- In this case, the conversion between 
angular apertures and transversal distances accounts for 
the redshift distribution of the GG sources. Values for 
gas overdensity are obtained from tt(0ap) estimates via 


^g(^Ap) 


{tt{9ap,z))^ 

(tTThe(^)T corr /(i±^))/ 


( 5 ) 


We find again good agreement with the dark matter pre¬ 
diction, finding no hints for feedback that is supposed to 
deplete gas from the inner halo regions. Similar results 
are obtained from Pksz(?') observations. 

Discussion. Observations of metallic lines in the X-ray 
in the direction of high energy sources provide informa¬ 
tion about the gas density in the so called circumgalactic 
medium of the Milky Way [1, 0, S [m • While some au¬ 
thors claim to have found evidence for all the baryons 
expected in our halo some other authors seem to 
find only between 10 to 50 % of the expected amount of 
baryons 0,S 1^. Both the thermal Sunyaev-Zeldovich 
[hereafter tSZ,[22| and the kSZ effects provide alternative 
approaches to detect ionised gas, different to the window 
available from X-ray observations. The tSZ and kSZ ef¬ 
fects provide statistical measurements of gas fraction in 
the visible halo population and do not restrict to the 
Milky Way host halo. In it is claimed that there is no 
apparent evidence for feedback effects in the tSZ lumi¬ 
nosity of GGs, in contradiction with X-ray observations. 
Those analyses have been recently revisited by 0 after 
including explicitly the impact of feedback in the filters 
extracting the tSZ signal. They conclude that Planck 
tSZ observations up to 5 R^qq are compatible with AGN- 
induced feedback effects, and incompatible with the no¬ 
feedback hypothesis (or self similarity in the tSZ luminos¬ 
ity - halo mass relation). Those results restrict however 
to gas collapsed in halos that is able to generate tSZ sig¬ 
nal 0. On the other hand, our kSZ study is blind to 
any assumed gas temperature profile in GG host halos, is 
not restricted to collapsed/virialised gas, and it provides 
lower limits to the amount of gas in GGs given the com¬ 
pensated structure of the AP filter. Figure [3] shows that 
our measurements are compatible with having detected 
all missing baryons in case these follow the dark mat¬ 
ter distribution. There exist uncertainties linked to the 
predictions for the gas peculiar velocities, but these are 
relatively small if we compare the predictions from the 
Millennium and the Gaussian simulations in the right 
panel of Fig. [TJ The impact of contaminants (such as 
dust or tSZ) seems to have been characterised and kept 
under control in [l^ , where it is also shown that a satel¬ 
lite fraction (of ^ 12 %) in our GGG should decrease the 
kSZ amplitude by less than 10 %. We thus conclude that 
the measured kSZ signal provides evidence for all the 
missing baryons predicted to be inside and around the 
GGs, which correspond to roughly half the total amount 
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of baryons present in the Universe at z ~ 0.12 under the 
SDSS DR7 angular footprint. In the future, a more de¬ 
tailed comparison with state-of-the-art hydrodynamical 
numerical simulations should shed more light on these 
results. 
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